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CHAFIER I 


Summary 


Recovery factor and friction factor data are presented here for air at 
a stagnation temperature of 110°F flowing through a cylindrical tube at 
Mach numbers ranging from 1.3 to 2.6 and at inlet diameter Reynolds 


numbers from 0.47 x 10° to 4.8 x 10°. 


Varying from a value of about 0.81 at a Mach number of 1.3 to about 
0.905 at a Mach number of 2,6, the recovery factor has been shown to be 
an almost linear function of only Mach number over the entire range of 
inlet diameter Reynolds numbers. Although the results obtained indicate 
three fiow regimes, there appears to be no dependency of the recovery 
factor on the diameter Reynolds number at inlet. In addition, there 
apparently is no simple relationship between recovery factors and Prandtl 
number involving fractional exponents. 

The mean apparent friction factors are in good agreement with the 


results of Keenan and Neumann!? and Ketchum! 9, From a plot of mean 





apparent friction factor as a function of length Reynolds number three 
flow regimes which are essentially laminar, transitional, and turbulent 
can be identified with the transitional regime apparently lying between 
length Reynolds numbers of 15 x 10° and 30x 10°. The friction factor 
attains a minimum value of about 0,0013 at the transition region and 
after transition decreases steadily with increasing length Reynolds 


number from a maximum value of about 00,0035, 


> aa 











Preliminary results were obtained with the recovery factor apparatus 
modified to reduce heat transfer to the test section; although neither ex- 
tensive nor conclusive, these data indicate that the recovery factors 
obtained will be from one-half to one percent lower than those determined 
with the unmodified apparatus and that recovery factor is independent of 


stagnation temperature, 
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CHAPTER Il 


Nomenclature 


Meaning 

cross-sectional area of test section (in*) 

cross-sectional area of nozzle throat (in~) 

specific heat at constant pressure (Btu/lb-F) 

nozzle flow coefficient (dimensionless) 

diameter of test section (in) 

diameter of nozzle throat (in) 

local apparent friction factor; friction factor which is the mean 
of the values for two adjacent stations in the test section 
(dimensionless) 

mean apparent friction factor; friction factor which is the 
mean of the values for the first station and another station in 
the test section (dimensionless) 
gravitational acceleration (ft/ sec“) 

mass velocity in test section (1b/ sec-ft7) 

mass velocity in nozzle throat (ib/sec -ft*) 

height of mercury in leg of McLeod gage (cm) 

height of mercury in capillary tube of McLeod gage (cm) 
»vroportionality constant (ft-lb/ Btu) 

specific heat ratio (dimensionless) 


thermal conductivity (Btu/hr-ft-F) 
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Symbol leaning 


longitudinal distance from exit plane of nozzle to a point in the 


t cated 


test section (in or ft) 
Lmax longitudinal distance from entrance to cylindrical tube (exit 
plane of nozzle) to a point where the Mach number is 1.0 (in or ft) 
M Mach number (dimensionless) 
static pressure; P,)., - barometer; Fg, - downstream stagnation 
tank; F sas- gasometer tank; P,)- upstream of orificejsP,- 
pressure differential across orifice; F,,.- upstream stagnation 
tank; P.,- wall of test section; Pyac- evacuated chamber 
surrounding test section (psi or psia) 
Prandtl number (dimensionless) 
gas constant (ft-lb/1lb-F) 
recovery factor (dimensionless) 
Reynolds number; Rep- based on test section diameter; 
Re; - based on distance from exit plane of nozzle; Re,- based 
on distance from nozzle throat (dimensionless); Rep*- based 
on nozzle throat diameter 
rnean Reynolds number; value which is a mean for two adjacent 
stations of test section (dimensionless) 
radius of station temperature thermocouple (in) 


radial distance between differential temperature thermocouples 


(in) 














Symbol Meaning 


et 





r radius of test section wall (in) 


w 
‘ temperature (°F); taw- adiabatic wall temperature; tg,- down- 
stream stagnation temperature; toas~ temperature in gasometer; 


t..- mean stream temperature; t - room temperature; 


m room 


t_- temperature at station thermocouple; St, - temperature 
S 
differential between differential temperature thermocouples; 


tug > Upstream stagnation temperature; 


u 

ge temperature (°Fabs); same subscript as for t 

w mass rate of flow (1b/sec); w,- mass rate of flow through 
orifice; Ws- isentropic mass rate of flow; w*- mass rate of 


flow through nozzle; w ~ mass rate of flow into gasometer 


gas 
x longitudinal distance from nozzle throat to a point in the test 
section (in) 
LD viscosity (lb/sec -ft); pam ~ viscosity evaluated at mean stream 


temperature; ,~s*- viscosity at nozzle throat 


DP dimensionless group, 4flmax 
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Introduction 



















The experimental work reported in this thesis is part of Project 
DIC 6418 sponsored by the Office of Naval Research to determine recovery 
factors, friction factors and heat transfer coefficients for air flowing at 
supersonic velocities. More extensive knowledge of the heat transfer 
coefficient is becoming increasingly important in the design of aircraft, 
aircraft coolers, missiles, compressor and turbine blading. However, the 
heat transfer coefficients as defined for air flowing at supersonic speeds 
cannot be determined without a prior knowledge of Mover factors, data 
for which is currently insufficient and inadequate. It is the purpose of 
this thesis to determine the magnitude of recovery factors and the 
parameters upon which the recovery factor depend in the flow of air at 
supersonic speeds through a cylindrical tube. In addition, data on closely 
Peiied friction factors are to be obtained incidental to the determination 


of recovery factors, 


Recovery Factor 
The coefficient of heat transfer for compressible flow is defined as 
the rate of heat flow per unit area per degree temperature difference 


between wall temperature and the adiabatic wall temperature: 


_ _ VA (1) 


lw - Taw 
The ternperature difference, Tw - Taw, iS now almost universally acceptec 


where any compressibility effects are encountered since its use yields 


i 
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values of h which are independent of the temperature difference at high 
subsonic and supersonic speeds. in earlier investigations, however, the 
heat transfer coefficients were often based on the temperature difference 
between Ty and Tm or To. 

Since Taw is defined as the temperature the wall would have if it 
were in equilibrium with the stream and thermally isolated from its 
Surroundings, i,e., an adiabatic wall, it cannot be measured when heat is 
being transferred by the wall, However, by determining the value of the 
recovery factor defined in terms of Taw and the stream properties found in 
an adiabatic test, Taw can be predicted under the conditions when heat is 
being transferred by the wall. 

The definition of the recovery factor used is the dimensionless 


quantity: 


Rt. (2) 


in which the stagnation temperature T, can be measured at a point in the 


stream where the velocity is zero; T can be measured in an adiabatic 


aw 
test; and the mean stream temperature T,, can be evaluated from a one- 


dimensional analysis of wall and stagnation pressure measurements, 


Since in steady, one dimensional flow of a perfect gas in a duct of 


constant cross-sectional area in the absence of heat exchange 











where V is the velocity of the mean stream at the point where T,, is 


evaluated, 





RS ae (4) 


jet ot eee (5) 
C 
2eJ p 


and T = T 


Equation (5) thus furnishes a method of predicting the value of Taw when 
heat is beins transferred by the wall, provided 2’ has been previously 
determined from adiabatic tests. 

From Equation (5), it can be seen that physically the recovery factor 
may be thought of as the proportion of the mean stream kinetic energy that 
is recovered at the wall as the air is retarded to zero valenae at the wall 
Surface. It is also seen that T,,, may be less or greater than T, depending 
upon whether R’ is greater or less than 1.0. 

In general, Taw will not equal T, since there is a velocity gradient in 
the fluid layers near the walls as the fluid is brought to rest which indicates 
also the presence of a temperature gradient. The temperature gradient in 
the fluid causes heat to flow from the layers near the wall to those further 
removed which are ata lower temperature. The velocity gradient, however, 
causes 5 transfer of work energy towards the wall in overcoming the viscous 
forces present in the boundary layer. The ratio of heat flow away from the 
wall to the work energy directed toward the wall determines whether T, y, 
will be greater or less than To. 

Since the flow of heat away from the wall is controlled by the thermal 


diffusivity % = k/p Cp and the flow of work energy toward the wall is 


~~ ee - 








controlled by the momentum diffusivity, or kinematic viscosity, V = A. it 
might be expected that the net effect is controlled by the ratio of X to / 
which is k/4 Cp, the Prandtl number. 

Pohlhausen“* in a study of adiabatic wall temperature of an incompress- 
ible fluid flowing over a non-conducting flat plate with a laminar boundary 
layer predicted that the recovery factor was a complex function of only the 
Prandtl number, 

Emmons and Brainerd! ® using a differential analyzer to solve the 
differential equations of compressible and incompressible flow over a flat 
plate with a laminar boundary layer found that the recovery factor was a 
function of the Prandtl number with some dependence upon the Mach number. 
For Prandtl numbers greater than 0.5 and Mach numbers less than 2.7, 


- ue 
they predicted that the relationship, R’ = (Pale would be valid to within 


two percent, 


19 20 predict that the recovery factor 


Eckert and Weise and Ackerman 


i 
is of the order of jee for a turbulent boundary layer. 


Recovery factors for subsonic flow of air in brass tubes were measured 
by Nicolai“! who reported values of 0.80 to 0.905 for Mach numbers ranging 
from 0,!] to 1.0, Czapek and Marcuse in similar measurements reported 
values of 0.&5 to 0.96, 

Ever’* measured surface temperatures on various types of probes 
placed in a supersonic stream and found that for Mach numbers between 
1.5 and 3.0, the recovery factors varied from 0.85 to 0.96 with increasing 


Mach numbers, 


a 














Ketchum, | ¥ who compilea the recovery factor data obtained by Junge 


lin previous theses on this 


and Mar golskee“ and Connors and Helfrich 
Project using apparatus similar to that used in this thesis, found that the 
recovery factor for air flowing ina lucite ney at Mach numbers from 
1.0 to 2.6 is a function of the Mach number only, providedthe flow is 
turbulent and indicated that the flow became turbulent at inlet diameter 
Reynolds numbers of 1.1 x 10°, Ketchum's values of recovery factor, 
0.88 at M = 2.53, 0.84 at M = 1.78 and 0.80 at M = 1.41 were about 3% 
lower than the recovery factors obtained by Eber. At inlet diameter 
Reynolds numbers less than ii x 10°, Ketchum found recovery factors 
that were somewhat erratic but which gave a mean line that for the 
higher Mach numbers agreed with the recovery factors obtained at the 
higher Reynolds numbers, but which at the intermediate Mach numbers 
gave recovery factors two to four percent higher than those obtained at 


the higher Reynolds numbers, At the low Mach numbers the two curves 


remained paraliel, 


Friction Factors 

The determination of friction factors in this thesis was undertaken 
to furnish data for the design of any future recovery factor or heat trans- 
fer apparatus and to supplement existing data on friction factors for air 
at supersonic velocities in a tube. 

The apparent friction factor or coefficient is an approximation of 


the true friction coefficient. The latter is defined in terms of the shear 











stress at the duct wall and is extremely difficult to measure. Previous 
studies indicate the desirability of using the apparent friction factor asa 
convenient form of evaluating results. It is defined in terms of the stream 
pressures, stagnation state and flow per unit area and may be evaluated 
from one-dimensional compressible flow functions for adiabatic flow at 
constant area from the Gas Tables.’ See Appendix B. 

Experimental results and literature on the subject of apparent friction 
factors in tubes is very limited. Rooceene data for supersonic flow 
indicated agreement at fully developed flow with the von Karman-Nikuradse 
relation between friction factor and Reynolds number for incompressible 
flow. Keenan and Neumann!” also obtain results that show that the friction: 
factor approaches that of incompressible flow at distances in excess of 50 
diameters from pipe inlet. In addition, they show that during the development 
of the boundary layer for a fixed value of L/D greater than 30, the mean 
apparent friction factor decreases as the Reynolds number increases, 

Ketchum, in his analysis of the friction factor data obtained by 
Helfrich and Connors, found that the apparent friction factors obtained in 
supersonic flow in textolite and lucite tubes were in good agreement with 


the results obtained by Keenan and Neumann. 











CHAPTER IV 


Inve stigation Procedure 


The investigation of recovery factors for air flowing adiabatically at 
supersonic velocities resolves itself into a determination of the values 


necessary to the calculation of recovery factors by the following formula: 


Ra = T _ WR See oe ln) 
a te ani REC oh 
In order to determine these values the work reported in this thesis 
was started under a plan, the basic features of which were: 
1. The use of dry air at various stagnation temperatures and over a 
wide range of stagnation pressures; 
2. The use of a converging-diverging nozzle to obtain supersonic 
velocities with a maximum Mach number of about 2.6; 
3. The passage of the air at supersonic speeds through a cylindrical 
tube with ee of about 41 and with facilities for measurement of 
wall temperatures and static pressures at regular intervals along 
its length; and 
4, The reduction to a minimum of heat transfer to the supersonic 
stream by providing thermal insulation in the forms of low con- 
ductivity plastic materials and of an evacuated chamber surrounding 
the lucite Re Raeical tube. 


In view of the results obtained by previous investigators whose work 


apparently led to the conclusion that the recovery factor is a function of the 





a 











stagnation temperature, it was considered desirable to conduct experi- 
mental runs at the following temperature levels: ae 8 65°, Elen and 130 cm 
In addition, in order to cover a reasonably wide range of Reynolds numbers, 
it was decided to vary the upstream stagnation pressure from about 8 psia 
to 80 psia, the upper limit having been imposed by the inability of a down- 
stream ejector to handle greater flow rates. In all runs critical flow 
through the nozzle was maintained by the use of a steam-driven air 

ejector which maintained very low pressures in the section of the apparatus 
downstream of the nozzle. 

In order to reduce the effect of condensation shock on the static 
pressure the specific humidity of the air entering the test apparatus was 
maintained at about 1.5 x 107°, It has been estimated that the presence 
of this amount of water vapor would cause a condensation shock which 
would affect the static pressures in the negligible amount of about 1 part 
in 100,000. !® 

A pressure tap was located in each of the eleven bosses spaced at 
two-inch intervals along the test section, and since the pressure gradient 
across the boundary layer can be considered negligible, the pressure 
recorded was essentially the mean stream static value. As described in 
Appendix B these pressure readings (P,,) were used to evaluate the Mach 
Number of the flow at each station, from which the T__/T,,, values could 


be determined. 


— — 








) Of the three thermocouples located in each boss one was used to 
measure the temperature at a point 0.089 inches radially outward from the 
tube wall and the other two were used to measure the temperature differ- 
ential between two points displaced 0.45 inches radially. From the 
differential value the radial temperature gradient (degrees per inch) was 
determined and the wall temperature was computed by extrapolation. In 
the absence of any heat transfer to the test section, the wall temperature 
thus measured would very closely approximate the adiabatic wall tempera- 
ture at that point. 

Since it was suspected that the previously obtained conclusion 
regarding the variation of recovery factor with the upstream stagnation 
temperature was caused by the fact that the tests were not actually 
adiabatic, it was deemed imperative that the heat transfer be rendered 
negligible. As a result the design and recommendations set forth in 
reference 4 were used in an effort to obtain the predicted temperature 
differential across the bosses of 0.001°F. 

Recovery factor runs were made at the following auegnatieei tempera- 
tures: BZ 65°, and 110°F; but preliminary calculations indicated that the 
results would be similar to those obtained by Helfrich and Connors, namely, 
that the recovery factor varied with stagnation temperature. It was, 
therefore, found necessary to disregard the twenty-one runs made at 32° 
and 65° and to devote further effort to the reduction of heat transfer. The 


} fourteen runs made at me: were retained and their results submitted for 





IV-4 


the reasons indicated in Chapter VI. 

Subsequent investigation indicated that the heat transfer to the test 
section could be substantially reduced by Se ciine the modifications to the 
apparatus as described in Appendix K. Preliminary calculations of 
recovery factors obtained with the modified apparatus indicate that the 
yalues for a stagnation temperature of 110°F will not be appreciably 
affected but that those at other temperatures will be modified in such a 
direction as to show that the recovery factor is independent of the stagna-~ 
tion temperature. A more complete discussion of these preliminary 
results can be found in Chapter VI; however, the final results will not be 
available until reported by the investigators who are continuing the work 


of this project. 
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CHAPTER V 


Description of Apparatus 


The apparatus used in this thesis is largely that originally designed by 


Klingensmith” and Wyant 


for the recovery factor test section which was designed mainly by Shoulberg 


I> and assembled by Junge and Margolskee@ except 


4 


and assembled by the authors with some changes, The major improvements 


of the present apparatus over that used by Junge and Margolskee and Connors 


and Helfrich! for determining recovery factors are: 


(1) 


(2) 


(5) 


(6) 


Use of lucite test section similar to that used by Connors and 
Helfrich but surrounded by a vacuum chamber to thermally in- 
sulate the test section and to reduce the time required to reach 
thermal equilibrium. 

Use of differential thermocouples to measure temperature 
gradients in the station bosses of the test section in order to 

more accurately determine adiabatic wall temperatures. 

Improved method of construction and installation of the test section 


~ 


thermocouples. : 
Use of a shorter nozzle to reduce the boundary layer effect upon 
determination of recovery factors. 

An improved method of joining the nozzle to the test section to 
reduce misalignments to a minimum. 

Alterations to the pressure measuring system to provide a method 


of measuring absolute pressures at wall of test section by means 


of DC7063 Silicone oil manometers, 











ers 


(7) Installation of a standard orifice and alterations to the downstream 
piping to provide sufficient straight lengths of uninterrupted piping 
for the orifice. 

(8) Installation of a source of dry air from the Low Temperature 
Laboratory. 

(9) Installation of a copper coil with a larger I.D. and a larger coil 
radius in the constant temperature bath to reduce the pressure 
drop in the coil. 

(10) Installation of several pressure gages and valves throughout the 
system to reduce the pressure losses and to aid in the control of 
the air and the testing of the system for leakage. 

(11) Installation of piping to permit calibration of nozzle and orifice 
by means of gasometer flow. (Piping is uncoupled at the apparatus 
when not in use and does not appear in any of the illustrations of 
this chapter.) 

A brief description of the components of the system appears in the 
following paragraphs with the exception of the orifice installation which is 
described in Appendix J. A recent modification to the test section, made 
after the data for this thesis was obtained, is described in Appendix K, 
Drawings and pictures of the apparatus, as it appeared prior to Run F-36 


and prior to the modification just mentioned, follow the description. 











Air Supply: 


Air is available from four Sources; 


Approximate 
Source Maximum Pressure Capacity 
(psi) (CFM) 
Ingerscll-Rand, steam driven, 2 stage 180 124.0 
Chicago ee ye lee 100 69.6 
Joy », motor a 200 250.0 
Low Temperature Laboratory 150 --- 


The last source is obtained from the excess dry air of the air liquefy- 
ing process in the Low Temperature Laboratory which uses the Joy air 
compressor for its air supply. Only air from the Ingersoll-Rand compressor 
or dry air from the Low Temperature Laboratory was used in the runs of 
this thesis. 

Air dehumidification: 

Refer to Figure !, Air from the Ingersoll-Rand compressor is passed 
through a coil immersed in an ice-water bath and a filter composed of glass 
wool and copper sponges to temove dirt and to condense oil and water vapor. 
The filter is equipped with a drain to vent the condensate to the atmosphere. 
The air then passes through a trap cooled by the ice-water bath and passes 
upward through a three-tube counter flow, regenerative cooler whose 
temperature approaches -60°F and in which the remaining water vapor is 
condensed out of the air in a solid phase, Leaving the top of the regenerator, 
the air enters a coil immersed in a dry ice-alcohol bath which cools the air 
several degrees, then passes through a filter of glass wool and copper 
Sponges where any solid phase is removed and returns downward through 


the regenerator. A thermocouple (T.) measures the temperature entering 
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the filter; and to prevent moisture condensation Shocks, air from the 
dehumidifier is not admitted to the test section when T, is higher than 
~40°F, During its downward path, the air removes heat from the upward 
flowing air and emerges from the regenerator as warm dry air. The de- 
humidification process is by~passed when air from the Low Temperature 
Laboratory is used since it has already been dehumidified in the air- 
liquefying process, 
Air Flow and Temperature Control: 

Air from the regenerator or from the Low Temperature Lab now 
passes throuvh the control valves which consist of a needle and a stop 
valve in parallel with a larger globe Stop valve. Fine control of the air 
flow is obtained by opening the large valve to the approximate upstream 
pressure desired and then using the needle valve to maintain the exact 
pressure, A by-pass valve is also provided just before the control valves 
to vent excess air to the atmosphere through a ile A shell-and-tube 
type Ross air heater, through which a mixture of steam and water is 
circulated, is used to heat the air to within approximately 5°F of the 
temperature of the constant temperature bath through which the air pro- 
ceeds next, This is a large water bath in which the air is heated to approxi- 
mately the desired upstream stagnation temperature as it passes througha 
copper coil (40° uncoiled length) immersed in the water, The water is heated 
by means of steam to the approximate temperature desired and the tempera- 
ture is then regulated by means of an electric heating coil in the bath with 


an automatic temperature regulator thermostatic control which operates to 








maintain the desired temperature within 0,2°F, Air is bubbled through the 
bath to stir the water and make the temperature uniform. After leaving the 
bath the air enters the upstream stagnation tank (3°' diameter) in which the 
air is slowed to a uniform low velocity by means of baffles and wire screen- 
ing. In the tank are nine thermocouples distributed in a plane perpendicular 


to the flow and equipped with radiation shields to prevent radiation to the 





tank walls. The thermocouples are used to measure the stagnation tempera - 
ture is) and the temperature distribution of the air in the stagnation tank. 
A pressure tap provides a means of measuring P,,.. 

Recovery Factor Apparatus: 

The air from the upstream stagnation tank enters the recovery factor 
apparatus through a convergent textolite reducer (Figure 3) which has an 
exit diameter of 125°", The reducer, upstream end plate and the nozzle are 
bolted together firmly with Allen head steel bolts. A stainless steel 
cylindrical ring (Figure 6) with a force fit in the upstream end plate connects 
the reducer to the nozzle, The ring is machined so that there is a near 
perfect transition with no discontinuities to interrupt the flow from the 
reducer to the nozzle. 

The convergent-divergent nozzle (Figure 4) is made of stainless steel, 
with an entrance diameter of 1.25", a throat diameter of 0.2685'’, and an 
exit diameter of 0.5015°’ and was designed for a Mach number of 2.6 and 
uniform parallel flow at the exit. A flange is provided at each end of the 
nozzle for securing it to the upstream end plate and reducer and for 


fastening the nozzle to the test section. The walls of the nozzle are as 


4 











v-& 


thin as possible but still provide the necessary strength. Two thermo- 
couple holes are provided to measure temperature distribution along the 
wall but the thermocouples were not installed, 

The nozzle, collar and lucite test section are bolted together firmly. 
The collar is provided to insure nearly perfect alignment of the nozzle 
exit and the test section entrance areas, Special machining and hand 
fitting was required in manufacturing the lucite collar to make the junction 
between the nozzle and test section as smooth as possible, Tne recovery 
factor apparatus was designed to be dismantled without the necessity of 
unfastening the nozzle and collar from the test section. Means of re- 
aligning the three parts is provided, however, should such an eventuality 
arise, 

The lucite test section (Figure 5) is similar to that used by Connors 
and Helfrich in regard to diameter (0.5015°"), length (21.020°"), and loca- 
tion of eleven test stations at 2°’ intervals along the tube. Lucite was 
used since it is a resin with a low thermal conductivity and reduces the 
effects of longitudinal heat transfer in the test section to a minimum. At 
each station, 1,.88'’-diameter bosses are provided to afford a means of 
anchoring the pressure tubing compression fittings to the lucite. In between 
the bosses the tube thickness is reduced to 1/8°’ to further restrict 
longitudinal heat transfer along the tube. In order to minimize the effects 
of turbulence on pressure measurement at further downstream taps 
the pressure taps are arranged around the tube so that only the first and 


eleventh pressure tap holes are in alignment. Three thermocouple holes 


SS 











are provided at each station. One thermocouple is used to determine the 





wall temperature while the other two, located at different radial distances 
from the tube wall, indicate the temperature pradient in the boss. The 
interior of the lucite tube was made as straight and smooth as possible 
and was givena high polish during manufacture, At the downstream end 
the test section was bolted to the downstream end plate which is also made 
of lucite. 

Surrounding the test section and bolted to the upstream and downstream 
end plates is a thick-walled outer tube (Figure 7) made of textolite for 
strength and low thermal conductivity. In the outer tube are two pipe con- 
nections to the vacuum system, used to evacuate the space between the 
outer tube and the test section, and a smaller pipe connection for the McLeod 
gage vacuum-measuring system. 

The upstream end plate made of lucite contains compression fittings for 
securing the pressure tubing and a means of leading thermocouple wires 
out of the apparatus. (See Figure 6). 

The recovery factor apparatus is sealed throughout with rubber gaskets 
and high vacuum grease at every joint to obtain as high a vacuum as possible. 


Glyptal and air-drying varnish are used externally to seal all joints and 





fittings. 


Downstream air path: 





The air leaves the recovery factor apparatus through a 4°’ I.D. sheet 
steel adapter fastened between the apparatus’ and the downstream stagnation 


tank. This tank like the upstream tank contains baffles and wire screens to 
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slow the air to a uniform low velocity, Nine thermocouples and a pressure 
tap are used to determine stagnation temperature (la) temperature dis- 
tribution, and Pq. Leaving the stagnation tank, the air enters a 2 1/2°° 1 
pipe leading to the orifice and from thence to the steam driven air ejector 
located in the Steam Laboratory. The air ejector reduces the pressure in 
the test section below atmospheric thus insuring supersonic flow through 
the greater part of the test section even at high rates of flow, 

Vacuum system: 

The vacuum system (Figure 1} used for evacuating the recovery factor 
apparatus consists of an all metal diffusion pump, a trap and a Cenco 
Megavac vacuum pump. The diffusion pump is a fractionating, water 
cooled, electrically heated vacuum pump using octoil and is capable of 
maintaining an ultimate vacuum of 1075 mm Hg. The Megavac pump 
requires a 1/2 hp motor and is capable of maintaining a vacuum of 
O.1 x io ° mm Hoe, A trap is provided between the two pumps to prevent 
accidenial flow of octoil into the Megavac pump. A flexible coupling is 
provided between the diffusion pump and the apparatus, which effectively 
reduces the vibrations transmitted from the Megavac pump to the apparatus. 
All piping joints are soft soldered and sealed with glyptal and air-drying 
varnish with the exception of two unions which are sealed with high vacuum 


-3 mm Reg. 


grease and glyptal. The best vacuum obtained was 1.81 x 10 
Pressure-measuring system: 


This system provides mercury and silicone oil (Dow-Corning 703) 


manometers as shown in Figure 10. Provided the pressure in the absolute 
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| manometer ranifold is less than 0,1 mm Hg, absolute pressures may be 
read directly from the silicone manometers which have one leg connected 
to the absolute manometer manifold which is evacuated by means of a 
Genco Hyvac Pump capable of maintaining 3 x 10-* mm Hg pressure. All 
silicone manometers are equipped with stop cocks to prevent oil from 
entering the pressure lines or the absolute manometer manifold. One 
manometer containing silicone oil is provided to compare the pressure 
in the absolute manometer manifold with that in the recovery factor 
apparatus. A McLeod gage calibrated by Junge and Margolskee is used 
to determine the pressure in the recovery factor apparatus, A Bourdon 


pressure gage is provided for measuring Py, greater than 43 psia. 


Silicone oil is used because of its low vapor pressure, low air solubility, 
and a specific gravity (1.09) near that of water. It also does not dissolve 
the grease used to keep the stop cocks air tight. 
Temperature-measuring system: 

The upstream and downstream stagnation tank thermocouples are 
those installed and calibrated by Junge and Margolskee, The thermo- 
couples installed in the test eee consist of eleven used to measure sta-~ 
tion wall temperatures and calibrated as in Appendix H and po ee 
used to measure the temperature gradient across each station boss. The 
thermocouples were made of #30 constantan and copper wire twisted 
together for a length of 1/8°’ and soft soldered. The junction was care- 
fully placed in the drilled hole of the station boss and the hole filled with 


high vacuum grease using a hypodermic needle. The 44 wires were led 





_— 
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out of the apparatus through 45° holes in the upstream end plate and thence 
between the end plate and the textolite outer tube, An effective seal was 
obtained by means of high vacuum grease, rubber gaskets, and externally 
applied glyptal and air-drying varnish. The cold junctions of the station 
wall thermocouples were placed in an ice water mixture contained in 
thermos bottles along with the cold junctions of the stagnation tank thermo- 
couples, The temperature reading was obtained by means of a Sensitive 
galvanometer and potentiometer in conjunction with a calibrated standard 
cell and a 6 volt storage battery. A switchboard was provided for selection 


of the desired thermocouple reading. 
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CHAPTER VI 


Discussion of Results 


General 

The results presented here were obtained from fourteen experimental 
runs, F-1 to F-9 and F-17 to F-21, using the apparatus described in 
Chapter V. Each run was made with the upstream stagnation temperature 
and pressure held constant while measurements were made at eleven sta- 
tions along the test section to obtain data necessary to the determination of 
recovery factors, friction factors, Mach numbers, and Reynolds numbers. 
All runs were conducted at an upstream stagnation temoverature of about 
110°F and at upstream pressures varying from 8 to 80 psia, the range of 
inlet diameter Reynolds numbers being 0.47 x 105 to 4.8 x 109. 

Runs F-17 to F-21 were carried out as duplicating runs to insure 
that the results obtained in Runs F-] to F-9 could be repeated. In general, 
there was excellent agreement between the data obtained from runs which 
duplicated each other, 

The results of the experimental work are presented in tabular form 
at the end of this chapter and are shown graphically in Figures VI-l to 
vI-10, Except for Figure VI-2 only data for points upstream ofa 
suspected shock were plotted. 

Figures ViI-! and VI~-2 show the distribution of wall pressures and 
calculated ediabatic wall temperatures along the test section for the 

} first nine oe From these plots can be seen the influence of snocks 


near the downstream end of the cylindrical tube as well as the effects of 


_ 
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heat transfer to the test section at its inlet and outlet ends. 

In the pressure-ratio plot, Run F-9, which was made with the lowest 
inlet diameter Reynolds number, shows a pressure ratio which increases 
along the length of the tube, a situation which is impossible in the one- 
dimensional theory used as a basis for the work of this thesis, Since 


16 


Similar results were reported by Kaye, Keenan, and McAdams’ ~ and 


10 it may be concluded that the discrepancy was caused by the 


Ketchum, 
limitations of the theory rather than by errors in measurement. 

In the temperature-ratio plot Run F-5 shows a distribution which is 
inconsistent with that of other similar runs, In checking through the 
original data used to calculate the results of this run, it was found that 
this was the first of five runs carried out in one day. As a result it is 
believed that temperature equilibrium had not actually been attained 


when readings were commenced and that, consequently, all results of 


Run F-5 should be disregarded. 


Recovery Factors 

The results which are of primary importance in the determination 
of recovery factors are those relating the recovery factor to the Mach 
number; they are shown in Figure VI-3, which includes the results of 
runs made over the entire range of inlet diameter Reynolds numbers from 
0.47 x 10° to 4.8 x 10°, Discounting the values for Run F-5, the deviation 
of any of the points from a mean line passed through them is less than 


one percent of the recovery factor, 
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In attempting to correlate the results presented here with those of 
previous investigators, a comparison with the data obtained by Ketchum 
was made. His data led to the apparent conclusion that the recovery 
factor depended on the inlet diameter Reynolds number as well as the 
Mach number; as a result he divided his results into two regimes, the 
line of demarcation being an inlet diameter Reynolds number of 
1.05 10°. Ketchum concluded that only in the region of higher inlet 
Reynolds could the recovery factor be regarded as a function of Mach 
number only. The results presented here do not indicate the necessity 
for such a division since the recovery factors determined seem to be 
independent of the Reynolds number. In addition, there is considerably 
less scatter in the region of low flow rates, 

At low values of Mach number the results indicated in Figure VI-3 
are in excellent agreement with those obtained by Ketchum in the regime 
of high rates of flow, but at the higher Mach numbers there is a signifi- 
cant deviation with Ketcitim's recovery factors being somewhat lower. 
The net result of this difference is that the mean of the values presented 
here is more nearly linear. 

In Figure VI-4 recovery factors are plotted against diameter 
Reynolds number for each run, The plot is interesting in that the lines 
connecting successive Stations for éach run are roughly parallel ata 
Slope of 0,47, This suggests a strong dependence of recovery factor 
on viscosity and, consequently, on the Prandtl number; however, it has 


been impossible to correlate the results indicated in Figure VI-3 with 
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Simple formules involving fractional exponents of the Prandtl number. 

The plot of recovery factor vs. length Reynolds number (Figure VI-5) 
shows two flow regimes which can best be characterized by the inlet 
diameter Reynolds number. If the results of Run F-9 are disregarded on 
the basis that the one-dimensional theory is inadequate for the conditions 
of the run, then one regime could be dened with flow for which the 


inlet Reynolds number is between 0,83 x 105 and 1.8 x 10° and for which 


the recovery factor could be expressed as a function of length Reynolds 
number, in the other regime for inlet Reynolds number greater than 
1.8 x 10° the recovery factor seems to be independent of the length 


Reynolds number and to be more dependent on the L/D ratio. 


Friction Factors 

Although local apparent friction factors were calculated along with 
the mean apparent values, the preliminary plots of the former showed 
no consistency and as a result, only the latter values have been plotted 
and will be considered in this discussion, The local apparent values, 
however, are included in the tables of results. 

The most significant results obtained are those indicated in 
Figure VI-6, which is a plot BE reietion factor as a function of length 
Reynolds ee. The plot is quite analogous to the corresponding 
relationship between friction factor and diameter Reynolds nurnber for 
incompressible flow. It clearly defines three flow regimes which appear 
to be essentially laminar, transitional, and turbulent. In this figure the 


region of transition lies roughly in the range of Reynolds numbers from 











15x 10° to 3ux 10°; this compares reasonably well with Ketchum’'s data 
which indicated that the limits of this regime were 7 x 10? to 25 x 10°, 

The friction factor reaches a minimum of about 0,0013 as it reaches 
the transition region and after transition decreases steadily with increasing 
Reynolds number from a maximum value of about 0.0035, 

In the plot of friction factor vs. L/D ratio (Figure VI-7) the results 
shown are in good agreement with the results of Keenan and Neumann! ? 
and Ketchum except for the values at the second station (L/D of about 5). 
In this report the frction factors for runs conducted at the higher inlet 
Reynolds numbers increase in going from station 2 to station 3 whereas 
in the other reports the friction factor decreases, Since the friction 
factor is dependent only on the Mach number, or the pressure ratio, it is 
not immediately evident why the difference exists. 

Figure Vi-8, which is plotted to a linear scale, again shows the 
presence of three flow regimes with a transition occurring at a diameter 
Reynolds number of about 10°. After transition the friction factor appears 
to level off at a value of about 0.003 or decrease eligneie from that value 
as the diameter Reynolds number increases, In addition, in the region 


beyond transition, the friction factor seems to be reasonably independent 


of the L/D ratio, 


Effect of Vacuum System 


Although no data was observed for the exclusive purpose of deter- 





mining the effect of the evacuated chamber on heat transfer, it is 


interesting to consider the apparent results obtained. Runs F-2 and F-18 


5 











were conducted at upstream stagnation pressures which differed by less 
than two percent and at upstream stagnation temperatures which differed 
by less than 0.1 percent, but F-2 hada FP... of 0.0035 mm Hg and F-18 of 
0.0091 mm Hg. In this case the ove factors of F-18 were roughly 
one-half percent higher than those of F-2. Unfortunately no other set of 
duplicating runs was made under such conditions as to afford a comparison 
as good as the above since in most cases the temperatures and pressures 
differec by greater amounts while the values of Pyac showed a smaller 


variation, 


Preliminary Results with Modified Apparatus 

After the recovery factor apparatus was modiiied to reduce heat 
transfer to the test section (See Appendix K), runs were continued in an 
effort to investigate the effect of stagnation temperature on recovery 
factor. The initial runs were made at stagnation temperatures of about 
110°F to check the validity of the results reported here. In Figure VI-10 
the results of two of these runs are compared with data obtained from 
the original apparatus for approximately the same diameter Reynolds 
number at Bee. It can be seen tone the recovery factor for the more 
recent runs is lower by about one percent at the higher Mach numbers 
and that at lower Mach numbers the difference is roughly one-half percent. 

The data obtained with the modified apparatus, consequently, more 
closely approximates that obtained by Ketchum despite the fact that his 


apparatus was more susceptible to heat transfer than that used for this 





report. However, since the differences are of sucha low order and since 


q 

















STATIC WALL PRESSURE = UPSTREAM STAGNATION PRESSURE 
RATIO vse RATIO OF DISTANCE FROM TUBE INLET TO 
TUBE DIAMETER 


Rey x 107 


Symbol _ Run_P : at inlet 
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Figure VJ-2 
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Run Ff «36 


Pus 
Tus 


Rep at inlet - 1,508 x 10° 


STN _ Pw/Pus 
1 .03970 
Z .04360 
3 .04589 
4 05287 
5 .06291 
6 .07138 
7 .08146 
8 <2935 

Run F -43 

Bis - 24.28 psia 

Poe — 265.02 OF abs 

ie - 00,0046 mm Hg 


- vac 


Rep at inlet - 1.457 x 10° 


STN 


l 
Zz 
3 
4 
5 
6 
4 
8 


- 24.91 psia 


Pyac - 0.011 mm Hg 


Py/Pus 


04490 





04088 


04745 
05383 
06314 
07084 
08259 
13151 


Taw/ Tus 


7247 
9314 
349 
9348 
.9368 
sJo9e 
9427 
Jolt 


R* 


8870 
tae 
8796 
8655 
855i 
6428 
63a 
15 he 


8864 
ise 
8 tam 
8652 
8528 
442 
8304 
of 507 
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Conclusions 


For air flowing in a tube at supersonic velocities the recovery factor 
is very nearly a linear function only of the Mach number in the range of 
Mach numbers from 1.3 to 2.6 and inlet diameter Reynolds numbers of 
0.47 to 4.8x 10°, The recovery factor varies from 0.81 to 0.905 at Mach 
numbers of 1.3 and 2.6 peated. 

Preliminary results with the apparatus modified to reduce heat 
transfer effects to a minimum indicate that the recovery means is inde- 
pendent of stagnation eae ae In addition the results indicate that 
the recovery factors measured with the modified apparatus are about 
0.5% lower than those reported in this ahi. 

Results obtained in this thesis show that it is possible to define 
three flow regimes in terms of length Reynolds number with a transition 
regime occurring between the Reynolds numbers of 15x 10° and 30 x 10°, 

In supersonic flow in a tube, the results show that variation of the 
mean apparent friction factors with length Reynolds number is analogous 
to the relation between friction factor and diameter Reynolds number in 
incompressible flow. The mean apparent friction factor reaches a 
minimum value of 0,0013 as it reaches the transition region and after 
the transition region it decreases steadily with increasing Reynoids 


numbers from a maximum value of about 0.0035, 
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CHAPTER VIII 


Recommendations 


Based on the results and experience obtained in preparing this thesis, 
the following recommendations are made: 

1, Further runs at stagnation temperatures other than 110°F should 
be made to determine whether the recovery factor is a function of the 
stagnation temperature, 

&. the range of Mach numbers covered should be extended by using 
a supersonic nozzle designed for a Mach number of about 4.0. 

3. In order to increase the range of flow rates which can be 
accommodated by the apparatus, the capacities of the ejector and re- 
generator should be increased. 

4. The actual effect on heat transfer of the evacuated chamber 
Surrounding the test section should be investigated to determine whether 
vacuums of the order obtained in this thesis are sufficient to render heat 
transfer necligible. 

5. in order to eliminate the effect of heat transfer to the air flowing 
between the stagnation tanks, it is suggested that a study be made of the 
possibilities of enclosing in an evacuated chamber the apparatus from the 
inlet of the upstream stagnation tank to the outlet of the downstream 
stagnation — it is further suggested that a means be provided for 
temperature control of the external surface of the test section. 

6, An attempt should be made to correlate the experimental data 


oy using Xeynolds numbers based ona viscosity which is evaluated at 


a 














Viale#é 





some temperature between the mean stream and the wall temperatures. 
?, For measuring upstream stagnation pressures greater than 
43 psig, a gage which is more accurate than the Bourdon gage now in 
use should be installed. 
8. The supersonic nozzle used in this investigation should be re- 


calibrated while it is not assembled in the recovery factor apparatus in 





order to better cover the range of expected nozzle throat Reynolds numbers. 
9. A filter should be installed in the pipe line through which passes 
the dry air supplied by the Low Temperature Laboratory. 

10, Safety traps should be installed in the low pressure legs of the 
absolute manometers in order to prevent the loss of silicone oi! through 
flow into the vacuum system manifold. 

11, The constant-temperature bath should be provided with a heater 
of greater power output since the one currently in use is incapable of 
maintaining a constant temperature of 110°F. 

12, A stop valve should be installed in the section of piping 
immediately upstream of the Cenco-Megavac pump to facilitate the 
detection of vacuum system leaks through the pump. 

13, The use of solid nylon as a material! for future test sections 
should be investigated. Mr. Chase, of the Cummings Machine Works, 
Boston, has shown that nylon in this form has a greater Pac haere ontsa 


‘than lucite, 














APPENDIX A 


Brief History of Project DIC 6418 


The work reported in this thesis was carried out as a part of 
Project DIC 6418, which is sponsored by the United States Navy through 
the Office of Naval Research. Since the Project has been in existence 
over a period of several years and since the results of previous investi- 
gators are directly related to those now being reported, a brief history 
is presented here, 

The following is quoted from reference 4, the publication of which 


immediately preceded this phase of the Project: 
‘‘The basic apparatus used to measure both recovery factors and 

heat transfer coefficients was designed in 1946 by Klingensmith (5). His 
work included the design of a supersonic (M = 2.5) nozzle, the design of 
a recovery factor test tube made of textolite to which the nozzle would be 
directly attached, and some calculations for the design of a heat transfer 
test section. The heat transfer test section was completely designed by 
Wyant (3) in 1947. 


‘Junge and Margolskee (6) began construction of the Klingensmith 
apparatus in February, 1947. The supersonic nozzle and textolite tube 
were machined at the Boston Navy Yard, The nozzle was calibrated by 
Ketchum (7) (who was in charge of experimental work between February 
1946 and November 1948) in early March, 1947, In early July, 1947, 
tests were made on the recovery factor apparatus. Prior to the assembly 
of the apparatus small burrs were visible in the test section at stations 
3and7, When test runs were made, the pressure distribution was not 
smooth, so the inside of the test tube was burnished with a glass rod to 
smooth the edges of the pressure taps, At the same time the test section 
was shortened by machining off two inches of the tube in order to give a 
greater length of supersonic flow, Further tests showed that station 5 
still gave an unsatisfactory pressure reading, so the hole was sealed with 
glyptal and a new pressure tap was drilled 1/8 in. further downstream, 
At the same time a brass collar was machined to fit between the nozzle 
and the test section. Additional trial runs were made during which it 
was observed that the pressure distribution was satisfactory near 
Station 5 when Pys was 21 psia, but unsatisfactory when Pys was 96 psia, 











The nozzle was then rotated 90° relative to the tube, but no appreciable 
change in pressure distribution was observed. Recovery factor runs 
JM1-9 were then made (6). 


‘*in October, 1947, Connors and Helfrich (4) took over operation of 
the apparatus and made runs T1 to T7 to verify the results of Junge and 
Margolskee (6). A lucite tube was machined and then installed. This tube 
had eleven stations and an extension which could be installed if desired. 
A diffuser was also provided, Since it was thought that the diffuser would 
increase the length of supersonic flow in the tube. The diffuser did not 
perform satisfactorily except at very low upstream Stagnation pressures 
(below 100 cm Hg absolute), At about the same time the manometer 
system was revised, 


‘*Runs LI to L.20 were then made (4), .with a range of stagnation 
temperatures between 33° F and 140° F, It appeared that the recovery 
factor might be a function of the upstream stagnation temperature, but 
careful investigation showed that the apparatus used by Connors and 
Helfrich was reliable only when tys was about 105° F, so no conclusions 
could be drawn about a Tys - R’ relationship. 

‘“In June, 1948, the heat transfer section was installed in place of 
the lucite recovery factor section by Larkin, Shoulberg, and Ketchum, 
At the same time the design of an improved recovery factor apparatus 
was begun.’’ 

During the period from June 1948 to May 1949 Shoulberg, Larkin, 
and England conducted an investigation of heat transfer coefficients for 
air flowing ina tube at supersonic velocities, Runs K~-1 to K-28 were 
made with the apparatus designed by Wyant. Runs 5-1A to S-3B were 
made after the apparatus had been modified to permit variation of the 
wall temperature without changing the air flow conditions. The super- 
sonic region investigated lay in the range of Mach numbers from 
1.0 to 2.5 and diameter Reynolds numbers from 30,000 to 460,000, 

In addition, Shoulberg, Larkin and England completed the design 


of an improved recovery factor apparatus which was subsequently 


modified and used for the experiments of this thesis. 














A= 


Working at the same time as the above, Ketchum, ina thesis for the 


doctorate degree, compiled the recovery factor and friction factor data 


obtained by Junge and Margolskee and Connors and Helfrich with the heat 


transfer coefficients obtained by Shoulberg and made a more thorough 


analysis of the data obtained. Ketchum also made a very complete review 


of the literature pertaining to this field. 


The work accomplished on the project from the period from June 1949 


to May 1950 by Welch, Volonte, Toong, and Young may be summarized 


chronologically as follows: 


l. 


2. 


EO. 





Design of Foelsch and Sauer supersonic nozzles. 

Revision of drawings and awarding of contract for the manufacture 
of the recovery factor apparatus. 

Making of subsonic heat transfer runs followed by dismantling of 
heat transfer apparatus. 

Assembly of new recovery factor apparatus without the test section 
thermocouples. 

Installation of pressure gages and valves in the air piping. 
Overhauling of regenerator and renewal of all filtering material. 
Design and installation of an orifice plate. 

Calibration of nozzle and orifice plate with the gasometer. 
Manufacture, calibration, and installation of test section 
thermocouples and wiring. 


Installation of larger copper coil in constant temperature bath. 























.and test section, 





14, Making of experimental runs F-36 to F-45 to determine the 


effects of the modification. 











APPENDIX B 


Calculation of Recovery Factors and Friction Factors 


Recovery Factors 

In general the method of calculation of recovery factors follows that of 
Helfrich and Connors.! By definition the recovery isi for a perfect gas 
is given by: 


Bee Taw ~ Tm a er Tus - Taw 
T 


Tus = a us a Tm 


Since T,,, is directly measurable and T can be closely estimated, it 


aw 
remains only to determine T,, in order to compute the recovery factor. 

For a gas flowing isentropically to a supersonic nozzle and adiabatically 
through the nozzle and duct of constant area, it can be shown that the product 
(P,,/P,) (A/A*) is the same as for the case of isentropic expansion, The 
pressure ratio, Mach number, and flow coefficient are related by the 


9 


expression 




















Since (A/A*) is fixed in this case, the Mach number can be determined when 

the flow coefficient and pressure ratio are known. Using tabulated values 

of (Py,/Py>) (A/A*) as a function of Mach number from the Gas Tables, ‘ 

a plot of (1/Cy,) (Py/Pus) val M was constructed and is shown in Figure C-1l. 
The pressure ratio at a given station can be evaluated from the observed 


data, and the flow coefficient can be determined from the nozzle calibration 

















plot of Cy, vs Rep* (Figure C-2) after Rep* has been computed from the 


following: 
_ (w/A*) D* 
Fas. 


Rep* where D* is a constant, 


From the isentropic tables for k = 1,4 

1+/ Tae 0.83333. 
Having computed T* the viscosity at the nozzle throat can be determined from 
the temperature-viscosity plot (Figure C-3). From the equation for maxi- 


mum isentropic mass rate of flow for a perfect gas through a nozzle: 


(w./A*) = 0.03012 _Pus __ 
us 


the isentropic flow rate is readily determined, In calculating the throat 
Reynolds number for the purpose of determining the flow coefficient, it is 
sufficient to use the isentropic rather than the actual value of flow rate. In 
an extreme case where a known flow coefficient of 0.9600 is neglected, the 
error in Repn* is Bee with a resulting error of 0.09% in estimated flow co- 
efficient, which error is well within the limits of precision with which the 
nozzle was calibrated. 

Once the Mach number has been determined from Figure C-l, the mean 


stream temperature can be evaluated from the following expression obtained 





from a one-dimensional analysis of adiabatic flow® 
ca |. 
7 K-1 We 
lis ® oO M 


The calculations are simplified by using the tabulation of the temperature 
ratio as a function of Mach number in the one-dimensional isentropic tables 


of the Gas Tables, 





er 











This method of recovery factor calculation differs from that of Helfrich 
and Connors in that: 

1, The flow coefficient has been included in the expression relating 
pressure ratio, area ratio, and Mach number; and 

2. In accordance with the assumption of a perfect gas the Joule~Thomson 
effect has been neglected in computing the stagnation temperature at each 
station, | 

The recovery factors for Run F-6 have been calculated by the methods of 


2 


Helfrich and Connors and of Junge and Margolskee;“ the results are presented 


for comparison in Figure Bel. 


Friction Factors 


From a one-dimensional analysis of compressible adiabatic flow ina 


tube of constant diameter it can be shown that? 
2 
4flmexs = 12oMi, Ket yj, (ke NM 2 p(M) 
D K M2 2K ae 


Since the @ function is tabulated as a function of Mach number for such flow 
in the Gas Tables, the calculation of friction factors is reduced to a deter- 


mination of the Mach number, which is evaluated incidental to the calculation 





of recovery factors, 
To compute the local apparent friction factor, i.e., the value which is a 
mean for two adjacent stations of the test section, it is necessary only to 


know the Mach numbers at the two stations. Then 


Afbmer - Afzlmoe 2 2, - Q 
O D 











Letting f be a value which is the mean of f, and fp 


Af laegeas real ae by, = G 
D 


But Lmax, ~ |-max2 is the distance between any two adjacent stations and 


is constant; therefore 


Ap ae = a0), = 0) 


D 
£ = b (p, - Dz) 
A al 





| 


One 


To compute the mean apparent friction factor, i.e., the value which is a 
mean between the first station and any other station, the above expression 


can be transformed to 


a 2 AUR. (>, ~- ®,) 


4 aly 


where a om is the distance between the stations involved. 
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APPENDIX C 


Calculation Data, Constants, and Coefficients 


Calculation Data 

Although the data used in the calculation of the recovery and friction 
factors was compiled in tabular form in order to insure consistency among 
the several persons making the calculations, they are appended here in the 
form of graphs. The necessary data include the following: 

l. Mvs (1/C,,) (Py,/P,.), the source of which is described in 
Appendix B and which is shown in Figure Cel. 

2. Cw vs Rep*, the source of which is described in Appendix G and 
which is shown in Figure C «2. 

aye vs T, which data is the same as that usec in references Il, 2, 4, 
and 10 and which is shown in Figure C-3, 

4. Fressure conversion data plotted in Figure C-4, 

5. Thermocouple calibration data, the basis for which is described in 
Appendix H and which is tabulated in Table C-l. 

6. T,/Tys vs M, which was obtained from the tabulation of T/ iene 
Table 30 of the Gas Tables (One-Dimensional isentropic Compressible-Flow 
Functions for k = 1.4). 


7. 4€L DC vs M, which was obtained from Table 42 of the Gas Tables 


mea af 


(Fanno Line - One-Dimensional Compressible-Flow Functions for Adiabatic 


Flow at Constant Area with Friction for k = 1.4). 














Calculation Constants 


The constants used in the calculation of recovery and friction factors 


are; 
A = 0.1978 in@ R= 53.345 ft-lb/1b-°F 
A* ~ 0,05662 in? r, - 0,340 in. 
D - 0.5018 in. Or, - 0.450 in. 
D* - 0.2685 in. ty - 0.251 in. 
g = 32.174 ft/sec¢ x,  - 2.076 in, 
k - 1.400 x2 - 4,076 in. 
| L,; - 0.810 in. . *2 = 6,076 in, 
Ly - 2.810 in. x4 = 8.076 in, 
Lg -4.810 in. x, ~ 10,076 in. 
Ly, - 06.810 in. my. = 12.076 an, 
Le + 8,810 in, x7 = 14,076 in, 
| Lg ~ 10,810 in, xg - 16,078 in, 
ia - 122810 in. xq ~ 18.076 in. 
Lg - 14.810 in, X)9 - 20.076 in. 
~ 16.810 in. 41, 7 22.076 in, 
- 18.810 in, Crabs = °F + 459.69° 


20.810 in. 








C) 
] 
te 


Calculation Coefficients 
Flow Rate 


The maximum isentropic mass rate of flow of a perfect gas is given by 














Z 


K 








Then Ws 2. 0.552 iene 
eo [Vas 
and Wie Ono SiG emis 
“Yus 
where w, is in Ib/sec, pys in psia, and Ty, in °Fabs. 


At the nozzle throat the isentropic flow rate per unit area is 


Ww 


: S 
G.* ome Pao 2543.4 Ws. 


where G.* has the units of Ib/sec-ft. At any station along the test section 


C.., W 


WwW Ss 
Gs ——p—— = 728.13 Cw we 


Diameter Reynolds Number 


The Reynolds number based on the diameter is given by 


Re 22 
D> = Sara 
an 
At the nozzle throat 
0.02337 G* 


and at any station along the test section 


0.04182 G 
an 


Ren _ 














where the G's have the units of lb/ sec -ft® and the viscosities which are 


determined at T* and T,,,. respectively, are in lb/ sec -ft. 


Adiabatic Wall Temperature 
The coefficient used in calculating the extrapolated adiabatic wall 


temperature is derived from the following: 


ts - taw : 7s ~ tw 
Ats . Ars ; 


For the constants of this report 


a 
tt 
‘N 


ts ~ taw 0.1978 Ats . 











Table C-} 
ihermocouple Calibration Data* * 





* * See Appendix H. 





a 


Mv tus&ty t & Ate 
-1.0 ~16.415°F 
~0.9 -11.415 
-0,8 - 6,451 
~0.7 - 1,523 
-0.6 3.58°F 3.369 
-0.5 8.40 8.225 
-0.4 13.18 13.045 
-0.3 17.94 17.828 
-0.2 22.65 22.576 
-0,1 27.34 27.306 

0 32.00 32.000 

0.] 36.64 36.648 
0.2 41,25 41,260 
0.3 45,83 45.854 
0.4 50.40 50.431 
0.5 54.94 54.97] 

0.6 59.45 59.493 
0.7 63.94 63.997 

0.8 68.41 68.465 

0.9 72.86 72.915 

1.0 77.28 77,347 

led 81,67 81.744 

lee 86.03 86.104 

1.3 90.37 90.446 

1.4 94.68 94.770 

1.5 98.99 99.076 

1.6 loseze 103,364 

7 107.54 107.635 

1.8 Ory, 111,869 

1.9 115.98 116.085 
2.0 120217 120,283 

zal 124.35 124.463 

22 128.51 128.625 

2.3 132.65 132.769 

2.4 136.77 136.895 

2.5 140.90 141.003 

2.6 144.99 145.112 

a 149.06 149.184 

2.8 153.12 153.238 

2.9 157.16 157.274 
3.0 161.19 161.292 
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Figure C-l 
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Figure C-2 
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Figure 0-4 
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Sample Calculation 


a 


Following is a tabulation of the operations involved in the calculation 
of recovery and friction factors, A sample data and calculation sheet for 
Run F -6 is included at the end of this appendix. The constants and coef- 


ficients used in the calculations are evaluated in Appendix C. 


Column Item Dimensions Operation 
(1) ay 8B none L/D values are constant; the 


values for stations 1-11 are, 
respectively, 1.614, 5.600, 9.586, 
Poul (6000, 21.092, 2.0 aes 
Z9.oe, 93.4909, 57.485, 41.471 
(2) x/D none x/D values are constant; the 
values for stations 1-1! are, 
respectively, 4.137, 8.123, 
12,108, 16,094, 20.080, 24.065, 


28.051, 32.037, 36.022, +40,.008, 


aoe 78 
(3) Troom oF Observed data. 
(4) Syne boreal ss Observed data (corrected to 32°F), 
(5) iS oe psi Multiply (4) by 0.49119. 
(6) Ee cm. Hg Observed data (manometer left 


column). 








Column [tem Dimensions _ a. Operation 
(7) Eas cm. Hg Observed data (manometer right 
column). 
(8) e " (7) ~ (6). 
(9) i psig Multiply (8) by conversion factor 


for mercury evaluated at 


temperature (3); Figure C-4. 


(10) 2 psia (9) + (5). 

(11) ey cm. Hg Observed data (manometer left 
column). 

(12) “ ms Observed data (manometer right 
column). 

(13) ‘ " (igen iy 

(14) = psig Multiply (13) by conversion 


factor for mercury evaluated at 


temperature (3); Figure C-4, 


(15) 2 psia (5) - (14). 
(16) a cm. Silicone Observed data (manometer left 
| column). 
(17) si a Observed data (manometer right 


column). 


(18) * ae (17) ~ (16). 











Column 


(19) 


(20) 


(21) 


(22) 
(23) 


(24) 


(25) 
(26) 
(27) 
(28) 


(29) 


(30) 


(31) 


Item 


Py 


us 


(1/Gy,) (Pw/Pus) 


Dimensions 


psia 


OF 


°F abs 
°F abs 


lb/ sec-ft x 1 o/ 


(°Fabs)@ 

lb/sec 

lb/sec-ft¢ 
none 


none 


none 


none 


Operation 
Multiply (18) by conversion 
factor for silicone oil 
evaluated at temperature (3); 
Figure C-4, 
Observed data; three values 
measured at different times. 
Convert average of three values 
in (20) using thermocouple 
calibration data, 
(21) + 459.69. 
0.83333 x (22). 
Find value of viscosity at 
temperature (23); Figure C-3. 
JV (22) 
0.03012 x (10)/(25). 
2543.4 x (26). 
0.02337 x (27)/(24). 
Flow coefficient evaluated at 
(28) from Figure C-2. 
(15)/(10) for mercury manometers. 
(19)/(10) for silicone oil 


manometers. 


(30)/(29). 








Column 
(32) 


(33) 


(34) 
(35) 


(36) 


(37) 
(38) 


(39) 


(40) 
(41) 
(42) 
(43) 
(44) 
(45) 
(46) 
(47) 
(48) 


(49) 


(50) 


[Item 


Tm/Tus 


ok 5 ca 
R’ 


ae Uae 


Dimensions 


none 


none 


OFabs 


oF 


si5h 
ots 
or 
oF 

OF abs 
none 
none 


none 


lb/sec-ft¢ 


lb/sec-ft x 10° 


x 10-9 


Operation 
Evaluated at (31) from Figure C-l. 
Evaluated at (32) from Table 30 
of Gas Tables. 
(33) x (22). 
Observed data, 
Convert (35) using thermocouple 
calibration data, 
Observed data. 
(35) + (37). 
Convert (38) using thermocouple 
calibration data. 
(39) = (36). 


0.1978 x (40). 


iG) Scie 
(21) - (42). 
(22) - (34). 
(43)/(44). 


] - (45). 


(42) + 459.69 
(22) 


728.13 x (26) x (29). 
Find value of viscosity at 
temperature (34); Figure C-3., 


0.04182 x (48)/(49). 








Column 


(51) 


(52) 


(53) 


(54) 


(55) 


(56) 


(57) 


(58) 


(59) 


Item 


Rep 


Re 


Re 


Re 


By - py 


D/4 aL, 





Dimensions 


x 10-2 


x 107° 


x 107? 


-5 


none 


none 


x 102 


none 


eS) 


rae te 


x 10 


cm 


em 


Operation 
Average of values for adjacent 
stations in (50). 
(1) x (50). 
Average of values for adjacent 
stations in (52). 
(2) x (50). 
Evaluated at (32) from Table 42 
of Gas Tables. 
In (55) subtract value for 
station 2 from value for 
station 1, station 3 from 
Station 2, etc. 
0.062725 x (56), 
in (55) subtract value for each 
station from value for station 1. 
Values are constant; for stations 
l-il they are, respectively, 
infinity, 6.2725, 3.1363, 2,090e 
1.5681, 1.2545, 1.0454, 0.8961, 
6.7841, 0.6969, 0.6273. 
(58) x (59). 
Observed data. 


Observed data. 
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APPENDIX © 


Design of Nozzle 


In the design of the improved recovery factor apparatus by Shoulberg,* 


the nozzle contour designed by Klingensmith? and used in the original 
recovery factor apparatus was retained with changes only in the dimension 
reference plane and in the mounting arrangements. 

In an effort to obtain a nozzle of shorter length than the Klingensmith 
nozzle and thus reduce the boundary layer effect upon flow through the noz- 
zle, two nozzles were designed in June-July 1949 by G. J. Van Wylen, 
research assistant to the Project, and R. V. Welch using methods devised 
by Kuno Foelsch!4 and R. Sauer!? for the design of axially symmetrical 
supersonic nozzles. 

Following the method of Foelsch, the complete nozzle contour from 
throat to exit was designed for uniform parallel flow at an exit diameter of 
0.5015’ and an exit Mach number of 2.8 with isentropic flow. It was esti- 
mated that the boundary layer effect in this nozzle would reduce the 
isentropic exit Mach number to an actual value of 2.6. After manufacture 
and installation of this nozzle in the improved recovery factory apparatus 
the Mach numbers measured at Station #1, whichis 0.81 inches from the 
nozzle exit plane, range from 2.68 at high to 2.33 at low rates of flow 
indicating that the average actual exit Mach number is slightly less than the 


desired value of 2.6. 








The Foelsch method linearizes the characteristic equations of axially 
symmetrical flow by comparing the conditions of flow in a cone with those in 
a nozzle. Mathematical expressions for the nozzle contour can thus be 
obtained. The design begins with an assumed source type flow OA (See 
Fig. E-1) which is altered in region AE to give parallel exit flow. The 
straight conical section AD assures conical source flow on entering 
region AE. The throat is designed by using the equation of continuity to 
determine the throat diameter and constructing TD as a circle tangent to OA 
at point D and parallel to the axis at throat T. Circular arc TC is a con~ 
tinuation of circular arc TD, and circular arc CI is drawn to produce an 
inlet nozzle diameter of 1.250 inches. 

The Sauer method is a part graphical, part analytical, solution of the 
characteristic equations of axially symmetrical flow assuming that the 
flow is isentropic, supersonic and irrotational; it uses successive approxi- 
mations which converge rapidly. Using the Foelsch contour from points T 
to A as a starting point, the exit portion of the contour was designed by the 
Sauer method for an isentropic exit Mach number of 2.8, The contour was 
then adjusted to obtain uniform parallel flow at an exit diameter of 0.5015’’. 
A more complete explanation of the Sauer method follows in Appendix F. 

Figure E=2 compares the nozzle contours obtained by the Foelsch 
and Sauer methods and also shows the Klingensmith nozzle used in the 
original recovery factor apparatus. The Sauer contour is about 10% 


shorter than the Foelsch, but the latter has smaller cross-sectional areas 








between the throat and exit planes. Since both contours were obtained by 
assuming isentropic flow in the nozzle, it was felt that the slightly longer, 
but smaller area, Foelsch contour would have less boundary layer effect 
than the shorter, but larger area, — contour. For this reason, 
Professor A. H. Shapiro, Department of Mechanical Engineering, MIT, 
recommended manufacture of the Foelsch nozzle. 

The nozzle was manufactured of stainless steel by the Cummings . 
Machine Works, Boston, Massachusetts, to tolerances of + 9005 on inside 
diameters and of £ .0001 on inside concentricity. 

It is interesting to note that the design of the nozzle contour by the 
Foelsch mathematical computation method required about 15 man-hours 


while that of the Sauer graphical-analytical method required in excess of 


160 man-hours. 
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ARKPENDIX F 


Description of Sauer Method of Nozzle Design 


The purpose of this appendix is to outline the Sauer method described in 
reference 13 which was used in the design of an axially-symmetric, super- 
sonic nozzle with uniform parallel flow at the exit and with an isentropic exit 
Mach number of 2.8. The reference does not fully explain the application of 
this method to this type of nozzle to those unfamiliar with the theory of 
characteristics. Therefore, for the future reference of those who might wish 
to design a similar nozzle, a brief description of the theory of the Sauer 
Method and a more detailed description of its application to this type of 
nozzle is included in this appendix. The nozzle exit contour found by this 
method appears in Appendix E., 

List of symbols: 
PD - potential function 
x ~ longitudinal cylindrical coordinate 


radial cylindrical coordinate 


ae 
t 


velocity of sound 


OQ 
0 


A 
‘ 


Mach angle 


oblique-angular coordinate system tangent to Mach lines at 
each point 


«6 
aS 


Wer - component of velocity vector parallel to F axis = = 
We = 29 99 99 99 2 aaa 2 = 2p 
ay 
Vv. 
ax 


= 


first approximation of Mach net 








M,; - second approximation of Mach's net 


Gy ~ first approximation of velocity net 


M 


second approximation of velocity net 


angle which velocity vector makes with axis in velocity net 


Part A-Theory: 

Flow is assumed to be steady, isentropic, supersonic, irrotational, and 
axially-symmetric, 

From the definition of the potential function d » equation of continuity, 
the Euler Equation, and the definition of sound velocity, there results the 


potential equation: 





Z 
do a sor 2 
so (Ge), 3 Gr) ]_2 3% 26 96) 1 265 
9 Xx? Cr an c2 c2 §axdy YX ar Lf ar Eq. (1) 
Assuming at each point in the physical plane a new coordinate system 
(E, n) tangent to the Mach lines at each point, the potential equation is 


transformed to: 


rp _ Sin 23> 
dB IN r ons | Eq. (2) 


If ina given region, we consider that this new coordinate system remains 


fixed while a given point P undergoes a small displacement, there results: 


a 
od 
aWep Wap _ Sind | oe 
mn a ¥ r 
For a finite difference, the atove may be written: 

sin of Eq. (4) 
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The significance of these quantities is shown in Figures F-1 and F-2. 
Consider a Mach net in the physical plane. A constant velocity is assumed 
in each Mach quadrilateral; hence this area will map in a point in the 
velocity plane. The center of each area is assumed to be at the intersection 
of the diagonals. A and my are the distances between the centers of two 
adjacent areas measured parallel to the ¥ and Y axis respectively, and are 
considered positive in the direction of flow. VY is the component perpendicular 
to the axis of symmetry of the mean velocity vector between the two adjacent 
areas, ris the distance from the axis to the center of the common side of 


two adjacent areas. 


Part B-Application: 

To illustrate the application of this method to the design of a nozzle, 
the steps in the design of a nozzle will be considered in some detail; values 
being taken from the design of this nozzle. 

I, Source flow ina conical nozzle. 

A. First approximation of flow in a conical nozzle: 

From the design by Foelsch method, the conical source flow 
was found to have a half cone angle of 11926° and the Mach number 


at point A was 1.875, which results in a Mach angle 


~| a 
%q_ = Sih a = p> (i 


For the first approximation in the physical plane, My Figure F -3, the 


Mach lines are drawn parallel to each other with an inclination of e a to 


the x-axis. Having this net, Ay ; oy and r can be determined. 











Points 1], 2, and 5 can be immediately located in the velocity plane, 
Figure F -4, since at each of these points the Mach number and the direction 
of flow is known. An initial approximation for point 3 is located by the 
intersection of characteristics through 1 and 2, which are inclined at +X A 
and - L respectively to the x-axis. From point 3 thus located, V for 2/3 


for 2/3 and A Wnp for 1/3 are calculated 


and 1/3 is determined, and A W : 


from Eq. 4, and point 3 is relocated. 

The initial approximations of subsequent points are determined by 
extending the characteristics of previously determined points and are then 
corrected in the manner explained for point 3. The points so located deter- 
mine the velocity plane G,. By measuring V and A from Gp by means of a 
subsidiary ellipse, Tene may be relocated once again to produce 
velocity plane Gy. 

B, Second approximation of flow ina conical nozzle: 

From Gy, the direction of the velocity vectors, 6, at each point 
is measured, and the Mach angle X is determined by using a sub- 
Sidiary ellipse. The direction of the characteristics at each area 
in the physical plane is given by 8 t£% . The new Mach net My] 1s 
constructed by drawing the new Mach lines with directions as 
determined by interpolating between the corrected directions. 
After locating the new centers of the quadrilaterals in the Mj] 


plane, they are located in the G,; plane to produce the G,' net. 
The points of the G,' net have the same relative position with 


respect to the points of the G, net, as do the centers of areas in 











of 


the M;, net have with respect to the center’areas in the M; net. 


To construct Grp AW and an are calculated from Mry and 


a 
Gy’. The corrected Mach directions found from G,; are drawn 
through the centers of areas of My. 

C. Third and higher approximations of flow in conical nozzle: 

F'or the third and higher approximations, the procedure for the 
second approximation is continued until for a given approximation 


there is negligible correction necessary to the Mach net in the 


physical plane, 


Il, To correct the conical nozzle for parallel outflow: 

The corrected nozzle results from joining a new Mach net M* to the 
Mach net M7, (or higher), so as to obtain parallel flow at exit with an 
isentropic Mach number My, = 2.8. The point of isentropic Mach number 
My. = 2.8 is located on the x-axis of the velocity plane of Gy,,;; and from 
corresponding characteristic lines is located in the physical plane of 
Mitt 

To construct M;*, draw characteristic a-a thru Mp. (See Figure F-5.) 
The centers of areas at fometn of My; and M; nets are to lie on line a-a, 
In My* the directions of the Mach lines are found by extending the character- 
istics of M,,;;. Upstream from a-a, the Mach net and velocity nets remain 
unchanged. Downstream from a~-@, the Mach and velocity nets are deter- 
mined by new boundary.conditions; namely, uniform parallel flow of 


Mr = 2.8 along a~b, 











r=-5 


A rough estimate of Go* is made remembering the new boundary 
conditions, 


calculated by means 


8 ° A 
Gy is cons(unet ae with values of ae p and W np 


of Mj* and Go*. Attention must be paid to the algebraic signs of db and a| - 
Construction must start at point Mr. 

With corrected directions of Mach lines in M,*, the second approxima- 
tions and higher are obtained in a manner Similar to that outlined for conical 
flow, 

When the correct Mach net is obtained, the contour of the nozzle is 
drawn. In each area through which the contour passes, the direction is 
M@ermined by the velocity vector for that area. The ratio of exit area to 


throat area (as determined by the Foelsch Method) should be that required 


to give the desired exit Mach number, 











axis of nozzle 
symmetry 


Fig. F-| Mach Net - Physical Plane 4 
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Fig. F-3 Physical Plane M™ 








Point 3 as determined by initia! ‘approximation 


Point 3 as calculated with values of V and 


determined from initial approximation 





Fig. F-4 Velocity Plane G, 
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Fig. F-3 Physical Plane Corrected to give 


Parallel Outflow 











APPERDiUte-G 


Nozzle Calibration 


Procedure 

Prior to calibrating the nozzle it was installed in the recovery factor 
apparatus in order that the calibration might yield flow coefficient data which 
included the influence of all components of the apparatus on the mass flow 
rate, Asa result the process was not strictly a *‘nozzle calibration."’ 

In order to determine the flow coefficient, i.e., the ratio of actual to 
isentropic mass rate of flow, twenty-three experimental runs were made 


using both the gasometer!! 


in the MIT Mechanical Fngineering Laboratory 
‘and a sharp~edged orifice to measure the actual flow rate. This appendix 
will deal only with the data obtained from the gasometer; the results obtained 
from the orifice, described in Appendix J, were not sufficiently consistent to 
be of any immediate use. 

Ten of the runs, NC-1 through NC-10, were made with apparatus set up 
as indicated in Figure G-l. In these runs high pressure air was heated in the 
constant temperature bath and passed through the nozzle, test section, and 
orifice to the sasometer, the rate of rise of which determined the actual mass 
flow rate. In the remaining eleven runs, NC-11 through NC -21, the apparatus 
was set up as indicated in Figure G~2; under these conditions air was passed 
from the gasometer which had previously been filled from an auxiliary air 
line, through the constant temperature bath, nozzle, test section, and orifice 


to a steam-driven ejector. In this case the actual mass flow rate was 


calculated from the rate of fall of the gasometer. 





The latter method more closely simulated the conditions under which 
the recovery factor runs were made since the ejector maintained low 
pressures in the section downstream of the nozzle. or this same reason 
the pressure ratio across the nozzle was always below critical and the con- 
ditions for maximum isentropic flow rate always existed. However, since 
the pressure of the air in the gasometer is limited by the number of weights 
which can safely be placed on the tank, the range of nozzle throat Reynolds 
numbers possible was very small and included only values which were 
generally below 10°. The range of Reynolds numbers was extended some~ 
what by uSing ice in the constant-temperature bath to reduce the temperature 
and viscosity of the air. 

In the former method the lack of a pressure tap in the nozzle makes it 
impossible to determine the isentropic mass rate unless the conditions for 
maximum isentropic rate exist. AS a result, it was necessary to experi- 
mentally determine the upstream stagnation pressure below which supersonic 
flow did not exist at the first station in order to determine the range of 
pressures for which the isentropic flow rate could be computed; it was found 
that the value was 61 psia. Since it is possible to simultaneously have 
subsonic flow at the first station (the first point in the test section after the 
nozzle throat at which there is a pressure tap) and critical flow through the 
_nozzle, there is a small range of upstream stagnation pressures for which 


the flow, though maximum isentropic in value, could not be computed. 





Calculations 

The calculation of the actual mass rate of flow is based on the assump- 
tion of a perfect gas. From the temperature and pressure of the air in the 
gasometer the density can be calculated and from the rate of rise or fall 
the volume rate of flow can be determined since the cross-sectional area of 
the gasometer is fixed at 19.79 square feet. The product of the volume rate 
and density yields the mass rate of flow. 

The maximum isentropic mass rate of flow and the nozzle throat 
Reynolds number are computed by the method outlined in Appendix B. 

A summary of the results of the calculations is tabulated in Table Gel 
and the plot of flow coefficient as a function of nozzle throat Reynolds number 


is shown in Figure G-3. 


Discussion of Results 

From Figure G-3 it can be seen that no calibration data was obtained 
for the Reynolds number range from 1.1 x 10° to 4.0 x 10°. Although it was 
possible to obtain flow at Reynolds numbers in this range, the physical 
limitations of the apparatus made it impossible to calculate the isentropic 
mass rate of flow. Asa result, it was necessary to take advantage of the 


data obtained by Ketchum! 9 in calibrating a similar nozzle. 


In carrying out the calculation of recovery factors the relationship 
between flow coefficient and nozzle throat Reynolds number determined by 
Ketchum was used. The straight line shown in Figure G-3 represents this 


function. 
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The results of the calibration carried out for this report are indicated 
by the points plotted in Figure G-3. It can be seen that the maximum deviation 
from the Ketchum line of any point is roughly 1.1% and that the mean deviation is 
considerably less, For this reason and in view of the fact that no better data 
was available at the time of the calculations, it is felt that the use of Ketchum's 


data is justifiable. 


Recommendation 
In order to determine the flow coefficient as a function of nozzle throat 
Reynolds number over the range for which no data was obtainable by means 
of the gasometer, the following procedures are possible: 
1. A pressure tap can be located in the nozzle such that the isentropic 
mass rate of flow.may be determined when the flow is not critical; or 
2. The nozzle can be removed from the recovery factor apparatus and 
calibrated independently of the apparatus. 
In view of the fact that it is undesirable to disturb the inner surface of 
the nozzle, it is recommended that the second procedure be carried out to 


calibrate the nozzle. 
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APPENDIX H 


Thermoc ouple Calibration 


Thermocouples used for the measurement of the upstream and down- 
stream stagnation temperatures ss and T a6) were those previously 
calibrated by Junge and Margolskee”. Results of this calibration are also 
found in Appendix C. 

Thermocouples for use in the lucite test section of the improved 
recovery factor apparatus to measure station temperatures and station 
radial temperature gradients were made of the same spools of #30 con- 
stantan and #30 copper wire as were used for the Tg and Tas thermo- 
couples. The manner in which the lucite test section thermocouples were 
made and installed in the apparatus follows that recommended by 
Shoulberg’. The station temperature thermocouples were calibrated in 
the MIT Heat Measurement Laboratory with the hot mation in a steam 
bath and the cold junction in at ice bath. The method of calibration follows: 

Symbols: Pj, ~ barometric pressure (corrected) 

T - temperature of boiling point of water °C 


mv__ = average potentiometer millivoit readings of the 
rp : : 
eleven stations while at steam temperature 


Trp - Bureau of Standards Conversion Table (for #30 
constantan and #30 copper wire thermocouples) 
temperature corresponding to mv,_, i.e., 

average temperature indicated by the eleven 
stations while at steam temperature 


mv_=- potentiometer millivolt reading at any tempera- 
ture between 0 and 100°C. 








ty - difference between Typ and T 


t. - temperature correction to be applied to t, 


t. - Bureau of Standards Conversion Table 
temperature corresponding to mv 


t, - true temperature at the station = t, ~.t. 

is 758.292 mm Hg 
T, = 100.0 +0.0369 (p,-760) -0.000020 (pp-760)2 = 99.937 °C 
MVry = 4,2884 mv 

- re) 
Trp = 100.27 °C 
t, = T_-T =0.33 °C 
L rp sp 


Assuming (a) that at 0 aC, the potentiometer reading (mv _) will be 
zero so that te = 0 and 
(b) that the temperature correction fe is a linear function 


of mv, then 


t i ween! s OF 
= oa. = = 1385 














mv; MV yp 4.2884 mv 
therefore t,=t,-—t. =t. - mv, (.1385) Op 
A table of values of t, vs mv was then computed by assuming 
values of mv, at 0.1 mv intervals, obtaining t, from the Bureau of Standards 
Conversion Table and then using the above equation to compute t,. (See 
Appendix C.) Linear interpolation between these values results in errors 
much smaller than the accuracy of the thermocouple measurements. 


At steam temperature the average deviation between the station 


millivolt reading and the averaged millivolt readings of the eleven 





stations (mv, ) ist 0.0025 mv or + 0.09 “= At temperatures between 


0 and 100°C, the average deviation will be proportionately less. 
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APPENDIX I 


Vacuum System Precautions 


The vacuum system itself is adequately described in Chapter V. The 
purpose of this enpendix is to give the en important precautions to be 
observed in the operation of the vacuum system in order to maintain a low 
vacuum of 0,003 to 0.005 mm Hp pressure, 

One of the major sources of potential trouble is the diffusion pump 
primarily because the pump is made of metal and its operating fluid cannot 
be observed nor can mal-operation of the pump quickly be detected by 
means of the vacuum measuring system. The condition and amount of 
octoil in the pump is not readily determinable since inspection requires 
breaking the seal on the not too accessible pump drain plug, removing the 
octoil and replacing the oil after measuring and inspecting. Inspection of 
the trap placed between the diffusion pump and the megavac pump to 
prevent carry over of octoil into the megavac pump is also difficult since 
its drain plug is quite inaccessible. Because of the danger of disturbing 
the Point eekla in nearby piping and also because about 5 hours of running 
time is required to remove all dissolved gases from any fresh octoil placed 
in the system, removal of the octoil from the diffusion pump for inspection 
is not recommended at intervals of less than two months unless some of 
the following precautions have been violated or the pump is definitely not 
operating satisfactorily. Because of the fractionating principle, however, 
one filling of octoil should last indefinitely if not subjected to abuse or 


action of chemical vapors. Therefore, the diffusion pump should be quite 





J Rad: 


dependable if the following precautions for the prevention of decomposition 

and contamination of the octoil are observed: 
1, Do not turn on the diffusion pump heater until the megavac pump 
has reduced the vacuum system pressure to less than 0.1 mm Hg 
because the hot oil vapors will decompose on exposure to air at 
pressures above 0.1 mm He. 
2. In order to condense the hot oil vapors before coming into contact 
with the air, continue circulating cooling water after turning off the 
pump heater for at least 15 minutes before opening the vacuum system 
to the atmosphere. a 
3. Open vacuum system to the atmosphere very slowly. 
4, When system is open to the atmosphere, the circulating water must 
not be turned on since water vapor will condense inside the pump and 
contaminate the oil, 
5. It is preferable to leave the vacuum system at pressures below 
atmospheric when system is not in use than to leave system open to 
the iii. 
6. Regulate the circulating water such that the cooling coil and entire 
surface of the diffusion pump beneath the coil are quite cold during 
operation. Water leaving the cooling coil should be slightly warm to 
the touch. 
7 When oetoil is to be removed for storage during long periods of 
pump idleness, it should be stored in glass stoppered, amber 


colored bottles. 
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8. Cleaning of the pump should be undertaken whenever the vacuum 
system is dismantled. Acetone should be used for this purpose, 
followed by air circulation, and then evacuation of the pump. Carbon 
tetrachloride must not be used since this solvent clings to the metal 
parts and is difficult to remove. 

The megavac pump requires no special care other than proper lubrica- 
tion of the motor, maintaining correct level of clean high vacuum oil in the 
pump, prevention of overheating, etc. The pump should give trouble-free 
and efficient operation for about 100,000 hours. 

The remaining causes of poor vacuum come under the heading of leaks 
through the numerous joints and seals in the vacuum system, Constant 
vigilance and Sereeon of glyptal and air-drying varnish to suspected 
joints is mandatory. It is preferable to apply the glyptal and varnish when 
the vacuum system is open to the atmosphere and time should be allowed for 
thorough drying before placing the system in operation, if possible. The 
covering of glyptal and varnish should be removed by means of glyptal 
thinner or benzene before dismantling a joint. When assembling any part 
of the vacuum system piping or test section, it is Secon to use every 
means possible, such as, cleaning of parts before assembly, use of high- 
vacuum grease on all rubber gaskets and screw threads, proper tightening 
of bolts and compression fittings, etc., to prevent future leakage. 

Some precautions to be observed in the vacuum measuring system are: 

1, Even though a safety trap is provided, care should be taken when 


comparing the pressure in the absolute pressure manometer menifold 





with that in the vacuum system to prevent silicone oil being forced into 





the recovery factor test section as a result of the pressure in the manom- 
eter manifold being much greater than that in the vacuum system. 

2. After measuring the pressure in the vacuum system by means of 

the McLeod gage, lower the mercury reservoir to its lower position to 
facilitate evacuation of that part of the apparatus, 

3. The mercury reservoir mentioned above should be lowered slowly 
to prevent an air bubble entering the mercury column from the mercury 


column air cushion and thus raise the pressure in the vacuum system 


appreciably, 











APPENDIX J 


Orifice 


Design and installation 

A standard orifice plate for measuring the downstream flow rate was 
designed and installed in the downstream Section of the recovery factor 
apparatus in order to provide a comparison with the calculated nozzle flow 
rate and to detect leaks in the apparatus, 

The orifice was designed in accordance with haemence 12 tor maximum 
possible flow rate occurring at the highest expected upstream stagnation 


pressure and the lowest expected upstream stagnation temperature, Thus: 


Wo = 9319 A* Pus = .5319 (,05662) 199 = 1355 Ib/sec 


Vj bus 492 
The area of the orifice to be installed in the 2 1/2°’-pipe with radius 
pressure taps was determined for the above flow rate from the following 
equation taken from Section VII of reference 12: 


SNE ak 29 Pies 
We — = Rae lo (1) 


in which the orifice upstream pressure (Pio) and the pressure differential 
GaP.) across the orifice are selected from a consideration of the manom- 
eters to be used for their measurement, From the area (A>) thus determined, 
the diameter (D>) was found to be about 1,5"". 

The orifice plate was then machined from 1/8”’ stainless steel plate 
with an orifice diameter of 1.498 inches and with a 45°-bevel on the down- 


stream face. The plate was installed between two 2 1/2-inch pipe flanges 








with a rubber gasket on each side of the plate, The length of straight piping 
provided upstream of the orifice was 22 pipe diameters while that down- 
stream was 10 pipe diameters, both of which are greater than the minimum 
required by Section IV of reference ld. The.upstream pressure tap was 
installed at a distance of 1 pipe diameter from the orifice plate and the 


downstream tap at a distance of 1/2 pipe diameter. 


Calibration 

The orifice was calibrated after installation, using air flow measured 
by the gasometer at the same time the nozzle was calibrated, The orifice 
flow coefficient was calculated from the following equation: 


Wgasometer 
0.668 ALLEY | Ro AP, 2) 
R Tos 


Results of the calibration are shown in Figure J-1, but it should be noted 


he = 


Po 
P10 
greater than 25 psia when T,,, 1s 


that this calibration is valid only for < 0.2 which corresponds roughly 





to upstream stagnation pressures (P,,,) 


near 110 OF, 


Calculation of Flow 

The method of calculating the orifice flow rate is adequately described 
on page 56, reference 12, There remains only the necessity of indicating 
how the values were measured, Referring to equation (1): P)q was obtained 
from the orifice upstream pressure tap reading converted to psia; AP, was 
obtained directly from the pressure differential manometer across the 


d tear 
orifice; IT} was assumed to be equal to Tgg; E, the are multiplier for 








thermal expansion of the primary element due to the pipe temperature, is 
equal to 1 at room temperature; and Y, the expansion factor, was calculated, 
knowing the ratio of orifice to pipe diameter and the ratio of 4P, to P). 
Results: 

As can be seen from Figure J-1, the orifice flow coefficients determined 
_from the gasometer flow vary as much as £ 5% from the mean value of the 
orifice coefficient. Assuming the gasometer flow to be correct, because it 
yielded nozzle flow coefficients that were quite close to values previously 
determined, the orifice flow determination is not accurate enough to permit 
the orifice to be used for the purpose for which it was designed. 

The orifice is limited to flow rates in which < 0.2 and thus for 
Pus greater than 25 psia when T,. 110 °F. This limitation can be improved 
by redesigning the orifice for the smaller flow rates. However, it was felt 
that re-design with a smaller orifice area would not improve the large 
variation of the orifice flow coefficient to within the t 0.5% necessary to 


make the orifice usable as a means of determining leaks or errors in the 


determination of nozzle flow rates. Therefore, re-design was not attempted. 
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APPENDIX K 


Modification to the Recovery Factor Apparatus 


After completion of runs F-1 to 35, the recovery factor test section was 
modified to reduce the heat conduction and radiation between the 1/§-inch O.D. 
copper tubing used for measuring the eg ran wall pressures, and the lucite 
test section, One end of the copper tubing is exposed to room temperature, 
while the other is fastened to the lucite test section. As discussed in 
Chapter VI, experiments with dry ice placed on the exposed ends of the copper 
tubing showed that an appreciable amount of heat was being transferred from 
the room to the test section via the copper tubing; or vice versa, depending 
upon the relation between the station wall temperature and room temperature, 

The modification consisted primarily of removing all copper tubing 
inside the recovery factor apparatus, leaving only a length of 1 1/4 inches at 
the upstream end plate and at the station boss compression fitting. The 
removed copper tubing was replaced with low conductivity, 1/8-inch plastic 
tubing. The ends of the plastic tubing were slipped over the short lengths 
of the remaining copper tubing and a tight joint was effected by means of 
small tube clamps and a coating of glyptal. This change effectively prevents 
all but small amounts of heat transfer by conduction and radiation between 
the room and the test section by the pressure tubing. 

To reduce the other possible radiation effects in the apparatus, such as 
that between the lucite test section and the textolite outer tube, aluminum 
foil was wrapped loosely around the test section, around each pressure tube, 


and then around the assembly of the pressure tubing and test section. Both 





end plates were covered with several layers of aluminum foil and a layer 
was loosely secured to the inside of the textolite outer tube. The recovery 
factor apparatus was wrapped with aluminum foil externally, as well. 

In order to determine the temperature surrounding the lucite test 
section, and also to determine the temperature gradient existing between the 
test section and the textolite outer tube, two thermocouples were installed 
inside the apparatus, The thermocouple wires were led out of the apparatus 
through four 459 holes drilled in the downstream end plate in a manner 
Similar to the thermocouple holes drilled in the upstream end plate, and 
thence through the joint between the downstream end plate and the textolite 
outer tube. A tight seal was effected by means of high vacuum grease 
applied to the wires, the end plate gasket and the joint between the end plate 
and outer tube, Externally, this joint was coated with glyptal dnd air-drying 
varnish, 

One of the above thermocouples, called T; , was located about 1/2 inch 
from the lucite test section between the bosses of stations #5 and 6, Two 
layers of aluminum foil separate the thermocouple from the lucite itself. 
The other thermocouple, Ty, was placed on the same horizontal plane and 
in the same longitudinal position as T; but was fastened by adhesive tape 
directly to the inside of the textolite outer tube, 

After the modification was completed, run F-36 was made at Tys = 
109.89°F and the recovery factors obtained were found to agree quite well 
with those previously obtained at similar upstream stagnation temperatures, 


This indicates that the modification does not appreciably change the data on 









recovery factors already obtained at T,,, near 110°F, This was the expected 










result since at this T,,, the wall temperature in the lucite test section is 
approximately equal to room temperature and thus only small heat transfer 
was possible via the copper tubing. At other stagnation temperatures, howeve 
the heat transfer before the copper tubing was removed was proportionately 
greater and this modification will result in much more accurate recovery 


factors, 
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APPEA Sok 


Effects of Measurement Errors 


The equipment used to measure observed quantities in the recovery 
factor and friction factor runs is the same as that used by previous investi- 
gators working on other phases of this project. Since the limitations on the 
precision of measurements have been discussed in references 2 and 4, it is 
the purpose of this appendix to present an analysis of the effects of the 
measurement errors on the calculated quantities, 

The recovery factor is dependent on the Mach number, which is a 
computed quantity, and the adiabatic wall temperature and upstream stagna- 
tion temperature, both of which are measured quantities, The Mach number 
in turn is a function of the local wall pressure, upstream stagnation pressure, 
and flow coefficient, 


Consider the equation (See Appendix B): 





aya). o/ E* 
Pus A* mM oe “ - 1 ms 


Using the technique of taking the logarithm of both sides of the equation and 


differentiating leads to: 


d (Pv / Pas) = daCw _ cl M { (K-1)M* dM 


7 2 a M  @0+simM) M 





for A/A* and k both constant. 
For small changes in variables the differentials in the above expression 


may be revlaced by delta-quantities to indicate small finite differences, To 





Lind 


study the effects of errors in flow coefficient only assume that there are no 


errors in the pressure measurements. Then 


OS: “s (1 ip K - | M°¢ \ aM 
Cx jee = wal = M 








Assuming k = 1.400 and rearranging: 





aM o ‘< + 2 M?) Be 2 
oM : ( + .4 M7?) (ry, 


The above equation expresses the interdependence of the percent error in 
calculated Mach number, percent error in flow coefficient, and Mach number; 
the function is plotted in Figure L-1, 

From the same original equation and ina similar manner the effects of 
errors in measurements of wall pressure and upstream stagnation pressure 


can be determined. The resulting error equations are: 











aM Ss (1) + eae AS es 

M (\ + 4M?) Fus 

Riv a .-(\ Se il 
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These functions are also plotted in Figure Le-l. 
To determine the effect of errors in Mach number on calculated recovery 


factors, the same technique is applied to the definition of recovery factor 


io EUR ee 
es (aes al) 





Since age oer is a function of Mach nurnber and since Day Ties is 


appreciably constant, the following can be derived: 


es a 1O ( = Tat ee) AM 


R Me |: = (\ eens rem » | M 





mor Tay/Ty,. = 0.94 





A WR a 0.6393 \ aM 
R Swat o.3\g M 


Figure L-2 shows the percent error in recovery factor as a function of 


percent error in Mach number for the above value of ee) ss 


From the definition of recovery factor it can be shown that 





‘B \2' Ei \ AN 
ee \ - (Vm / er) TV ae 


This function is plotted in Figure L-3 for various values of Mach number 
and for Taw/Tus constant at 0.94. 
As a result of this analysis the probable errors in calculated recovery 


factors due to the limitations on the precision of measurements are tabulated 


> 





below: 
Probable Maximum Maximum 
Measured Measurement Probable Probable 
Quantity Error Error (%) R’ Error 
M= 1 M = 2a 
Bis (Bourdon gage) 0.5 psi 1.3% 1 %o 0.1% 
Be (Hg manometers) 6.10 cm 0.2% 0.2% 0.05% 
P., (Hg manometers) 0.10 cm 2% 1.6% 0.3% 
fe (absolute manometers) 0.10 cm 0.5% 0.4% 0.05% 
Cy 0.01 1% 0.8% 0.1% 


Taw Miaa/ Tus = .94) 0,4°F 0.1% 0.8% 0.2% 
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